Transmutation of momentum into position in magnetic vortices 
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We show that transmutation of linear momentum into position may occur in a system of three 
magnetic vortices thanks to a direct link between topology and dynamics in a ferromagnet. This 
happens via exchange between the linear momentum of a vortex-antivortex pair and the position 
of a single vortex during a semi-elastic scattering process. Vortex polarity switching occurs in the 
case of inelastic collisions. 
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I. INTRODUCTION 

Magnetic vortices have been discussed traditionally in 
the context of two-dimensional ferromagnets [j], 0, 
They are magnetization configurations with a nontriv- 
ial topological structure and are thus characterized by a 
topological invariant. Interest in vortices has increased 
dramatically in recent years due to many related observa- 
tions in mesoscopic ferromagnetic elements, such as disc- 
shaped thin elements with submicron dimensions 0, 0, [f| 
where vortices are created spontaneously. The magneti- 
sation vector is tangent to the side boundary of the disc 
and points along the disc axis at the disc center. Simi- 
larly, an antivortex is a configuration where the magneti- 
sation vector winds around the center in a sense opposite 
to that of a vortex and has been observed in a specially 
designed magnetic element Q , and as a metastable state 
in a ring element 0]. A significant feature of a vortex 
or antivortex is that its energy increases logarithmically 
with the size of the element and would thus tend to in- 
finity on an infinite film. This fact may explain why 
experimental observation of a single vortex is relatively 
recent. 

It is also possible to construct nontrivial magnetic con- 
figurations in the form of vortex-antivortex (VA) pairs. 
Unlike single vortices, VA pairs are localised configu- 
rations and have finite energy even on an infinite film. 
Vortex-antivortex pairs have been predicted to be created 
by alternating external fields in an early study based on 
collective coordinates @ . They play a central role in dy- 
namical processes that lead to vortex polarity switching 
in ferromagnetic elements as reported in recent experi- 
ments [13, til El- 

A subtle feature of a vortex is its polarity, i.e., the di- 
rection of the magnetisation at the vortex center, usually 
called the vortex core. Therefore, in a VA pair the vortex 
and the antivortex may have either the same or opposite 
polarities. In the case of like polarities, a VA pair under- 
goes Kelvin motion in a direction perpendicular to the 
line connecting the vortex and the antivortex [13j |, with 
velocity that is inversely proportional to the distance be- 
tween the two vortex centers. Such a VA pair will be 
called a Kelvin pair. In the case of opposite polarities, 
a VA pair undergoes rotational motion around a fixed 



guiding center, with angular frequency that is inversely 
proportional to the square of the distance between the 
vortex and the antivortex Such a VA pair will be 

called a rotating pair. As it turns out, the three-vortex 
process that leads to vortex-polarity switching involves 
both types of VA pairs and may be thought of as a com- 
position of Kelvin and rotational motion. 



Thus our aim in this paper is to study in some de- 
tail a three- vortex process in which a Kelvin pair collides 
against a single vortex initially at rest. During collision 
a transient rotating VA pair forms and its dynamics is 
crucial for the final outcome of the scattering process. 
In short, when the speed of the incoming Kelvin pair 
is sufficiently low, the scattering process is mediated by 
the formation of a rotating VA pair but is eventually 
semi-elastic; the Kelvin pair reemerges after collision at 
a peculiar scattering angle while the initial single vortex 
comes to rest at a new position. This unusual behavior 
can be explained by the special nature of the conserva- 
tion laws of linear and angular momentum due to the 
underlying nontrivial topological structure. In contrast, 
for larger values of the speed of the incoming Kelvin pair, 
the intermediate rotating VA pair suffers a topologically 
forbidden annihilation (which is possible on a discrete 
lattice) and thus leads to a highly inelastic scattering 
process; the final product is a single vortex at rest with 
polarity opposite to the polarity of the original single 
vortex. 



Numerical simulations are performed on the basis of 
the two-dimensional (2D) Landau-Lifshitz (LL) equation 
briefly described in Sec. II. In Sec. Ill we present our main 
result for semi-elastic scattering at low velocities, which 
may be understood as a transmutation of momentum into 
position (and vice versa) due to the underlying topol- 
ogy. Vortex polarity switching occurs at higher velocities 
and is described in Sec. IV. Some concluding remarks 
are included in Sec. V. In the Appendix we present an 
analytical solution based on collective coordinates which 
provides a reliable approximation at sufficiently low ve- 
locities, in agreement with the numerical simulations of 
Sec. III. 
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II. MAGNETIC VORTICES 



topological vorticity 7 = 7(2;, y, t) defined from [19|, [2C 



A ferromagnet is characterized by the magnetization 
m = (mi, ?7i2, 777.3) measured in units of the constant sat- 
uration magnetization M s . Hence m is a vector field of 



unit length, m 



1, but is otherwise 



a nontrivial function of position and time m = m(r,t) 
that satisfies the rationalized Landau-Lifshitz (LL) equa- 
tion 



= m X f , f = Am 

dt 



qm 3 e 3l m 2 = 1. (1) 



Here distances are measured in units of the exchange 
length £ C x = y/ A/2ttM^, where A is the exchange con- 
stant, and the unit of time is To = 1/(47T7M S ) where 7 
is the gyromagnetic ratio. Typical values are £ cx ~ 5nm 
and To ~ lOps which set the scales for the phenomena de- 
scribed by Eq. <JTJ) . To complete the description of the LL 
equation we note that we consider ferromagnetic materi- 
als with uniaxial anisotropy. Then e 3 in Eq. fT]) is a unit 
vector along the symmetry axis and the dimensionlcss 
parameter q = K/2ttM 2 , where K is an anisotropy con- 
stant, measures the strength of anisotropy. In particular, 
q is taken to be positive throughout this paper, a choice 
that corresponds to easy-plane ferromagnets. Also note 
that we have neglected the demagnetizing field in Eq. (JTJ) 
which amounts, in some respects, to a simple additive 
renormalization of the anisotropy constant in the thin- 
film limit 



involved 
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To be sure, the thin-film limit is more 



17L [IH but some of the key issues can be 
discussed already within the simplified model defined by 
Eq. (JTJ. We adopt this limit in the remainder of this 
paper and thus consider a two-dimensional (2D) restric- 
tion of Eq. (1) while q is set equal to unity by a suitable 
rescaling of the space-time coordinates x,y and t. The 
effective field / in Eq. ([1]) may be derived from a varia- 
tional argument: 

/ = -Jj|; E = \ J [( Vm ) 2 + m l] dxd y ( 2 ) 

where E is the conserved energy functional. 

Static solutions are obtained by solving Eq. (1) with 
the time derivative set equal to zero or, equivalently, by 
finding stationary points of the energy functional. We 
can write the general axially symmetric vortex solution 
as 



mi + i 777,2 



a iK,(<t>-4> Q ) 



7773 



A cos ( 



(3) 



where p and <p are the usual cylindrical coordinates and 
9 = @ (p) is the vortex profile that can be found numeri- 
cally [1 31 ] - The vortex number is re = +1 for a vortex and 
re = — 1 for an antivortex. The integer A = ±1 defines 
the vortex polarity. The constant <po is an arbitrary an- 
gle that reflects the azimuthal symmetry of Eq. (1) and 
does not affect the vortex energy. 

A key quantity for describing both topological and dy- 
namical properties of the 2D LL equation is the local 



7 = 2 e "/3 i d rn x d a m) ■ m, 



(4) 



where the usual summation convention is invoked for the 
repeated indices a and (3, which take over two distinct 
values corresponding to the two spatial coordinates x and 
y, and e a p is the 2D antisymmetric tensor. In particular, 
one may consider the total topological vorticity T and the 
Pontryagin index or skyrmion number Af defined from 



7 dxdy, Af 



r 

An' 



(5) 



For the vortex given by Eq. © the total vorticity T and 
the skyrmion number Af are calculated from Eq. ([5]) to 
be 



F=-2ttkA, Af = 



1 



(6) 



We must thus distinguish four types of vortex states la- 
belled by (re, A) where the vortex number re = ±1 and 
the polarity A = ±1 may be taken in any combina- 
tion. In all cases a vortex state is characterized by half- 
integer skyrmion number Af, in contrast to magnetic bub- 
bles observed in easy-axis ferromagnets which carry in- 
teger skyrmion number. Nevertheless, a magnetic vortex 
shares with a magnetic bubble the fundamental property 
that it is spontaneously pinned around a fixed guiding 
center; i.e., it cannot move freely unless an external mag- 
netic field gradient is applied or other vortices are present 
in its vicinity. 

One may also consider a pair of vortices («i,Ai) and 
(«2, A2) with total skyrmion number that is an integer: 



Af = — ~~ (reiAi + re 2 A 2 ). 



(7) 



The energy of such a configuration is finite if we restrict 
attention to vortex-antivortex (VA) pairs («i = — re 2 ). 
For definitcness, we choose rei = 1 and «2 = — 1 and thus 
the skyrmion number 



Af 



1 



(Ai - A 2 ) 



(8) 



depends on the polarities Ai and A 2 . Now, a VA pair with 
equal polarities (Ai = A 2 = ±1) is topologically trivial 
(Af =0). In contrast, a VA pair with opposite polarities 
is topologically equivalent to a skyrmion (Af = 1, for 
Ai = —1 and A2 = 1) or an antiskyrmion (Af = — 1, for 
Ai = 1 and A2 = — 1). 

Unlike single vortices, the VA pairs sketched above 
cannot be static solutions of the LL equation. A topolog- 
ically trivial (Af — 0) VA pair undergoes Kelvin motion 
in a direction perpendicular to the line connecting the 
vortex and the antivortex [l3j], whereas a topologically 
nontrivial (Af = ±1) VA pair undergoes a rotating mo- 
tion around a fixed guiding center [l4| . The gross features 
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of these two types of motion of a VA pair had been an- 
ticipated by a treatment based on collective coordinates 
[9] which is valid only in the limit of widely separated 
vortices. In contrast, Refs. [13] and [14] provide a com- 
plete field theoretical calculation within the LL equation 
of true solitary waves in Kelvin and rotational motion, 
respectively. Such a calculation remains valid for any 
relative separation of the VA pair, even when the pair 
reduces to a solitary lump with no apparent topological 
features and is thus beyond the reach of collective coordi- 
nates. Both types of motion of a VA pair are important 
for understanding the three-vortex process analyzed in 
the continuation of this paper. 



III. VORTEX PAIR SCATTERING 

Recent experiments have shown that a vortex may 
switch its polarity under the influence of a very weak 
magnetic field of the order of a few mT [l(], EH- The 
key to this phenomenon is the appearance of a vortex- 
antivortex (VA) pair which is spontaneously created in 
the vicinity of a preexisting vortex. We are thus moti- 
vated to study in detail this particular system of three 
vortices. Here we do not address the question of how a 
VA pair is actually created. In fact, production of a VA 
pair might not be possible in the idealized model stud- 
ied in the present paper. Rather we concentrate on the 
second step of the process and explain in detail how the 
vortex interacts with the VA pair. 

It is obvious that the full dynamics of a three-vortex 
system is rather rich. In order to be able to analyse 
this system in terms of elementary processes we set up 
the initial configuration as depicted in the first panel 
of Fig. [T] Specifically, we assume that a single vortex 
(k, A)=(l, — 1)=C is initially at rest at some specified 
point which is taken to be the origin of the coordinate sys- 
tem. We further assume that a Kelvin pair (with JV = 0) 
consisting of a (1, 1)=A vortex and a (—1, 1)=B antivor- 
tex is somehow created in the neighborhood of the origi- 
nal vortex. Once created the AB pair will undergo Kelvin 
motion along the y axis and eventually collide with the 
single vortex C. 

One could use any reasonable VA pair ansatz for the 
pair AB. Here, we shall use the VA pairs calculated in 
Ref. [13j as exact steady state solutions of the LL equa- 
tion. These are viewed as solitary waves which propagate 
with a velocity in the range < v < 1 where the maxi- 
mum velocity v — 1 coincides with the low-wavenumber 
limit of the magnon group velocity expressed in ratio- 
nalized units. For small velocities v <§; 1 these solitary 
waves are composed of a vortex A and an antivortex B 
at a relative distance d ~ 1/v. 

The process was simulated by integrating numerically 
the corresponding initial-value problem in the LL equa- 
tion. Figure [1] provides an illustration with three charac- 
teristic snapshots in the case of a relatively slow Kelvin 
pair initially moving along the y axis with velocity v = 
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FIG. 1: Three snapshots for the collision of a VA pair in 
Kelvin motion (the AB pair), initially located at (0, —15) and 
propagating with velocity v = 0.1 along the y axis, against 
a target vortex C initially located at the origin. During col- 
lision, antivortex B rotates around vortex C before rejoining 
its original partner A to form a new VA pair that scatters 
off at an angle in the third quadrant. The target vortex C 
is shifted to a new position in the fourth quadrant thanks 
to transmutation of VA pair momentum into vortex position. 
Left column depicts the 2D projections of the magnetization 
(7711,7712) while the right column depicts the level contours 
of 7773 with solid lines corresponding to positive values and 
dashed lines to negative ones. 



0.1 for which the vortex and the antivortex are separated 
by a distance d « 1/v = 10. As the pair approaches, the 
original C=(l, —1) vortex teams up with the B=(— 1, 1) 
partner of the AB pair to form a new, topologically non- 
trivial (J\f =1) VA pair in quasi-rotational motion. In 
fact, B rotates almost a full circle around C before re- 
joining its original partner A. The new AB pair is again 
a topologically trivial (AT = 0) VA pair in Kelvin motion 
that moves away from the target vortex, having suffered 
a total scattering angle that is greater than tt/2 from 
its original direction. The scattering is slightly inelastic 
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in the sense that the outgoing AB pair moves out with 
greater velocity (v — 0.15). And, most remarkably, the 
target vortex C moves away from the origin and comes 
to rest at a new location in the fourth quadrant of the 
xy plane. 

The unexpected result in this numerical simulation is 
that the VA pair AB is scattered at an angle to the y 
axis which is the initial axis of propagation. In a naive 
approach to the problem one would think that the VA 
pair scattering and the corresponding change in its lin- 
ear momentum does not seem to be accompanied by a 
corresponding change in momentum of any other object 
taking part in the process. This would appear to contra- 
dict the conservation of linear momentum in this conser- 
vative system. 

Such a peculiar behavior can be explained by the un- 
usual nature of the linear momentum which allows for a 
transmutation between position and impulse in the case 
of topologically nontrivial systems, such as the three- 
vortex system considered here. The two components of 
linear momentum (impulse) P = (P x , P y ) can be written 
as moments of the topological vorticity 7 [l9l |20T | : 

Px = - j ' m dxdy, P y = J £7 dxdy, (9) 

while the angular momentum is given by 

L=^Jp 2 1 dxdy (10) 

where p 2 — x 2 + y 2 . 

We first consider the impulse of a single vortex like 
vortex C shown in the first panel of Fig. Q] which sits 
at the origin. The impulse is clearly zero due to the 
axial symmetry of the vortex. It should also be noted 
that P is actually a measure of position for an isolated 
vortex. Indeed, if the vortex is shifted to a new posi- 
tion {x, y) = (a, b), the corresponding impulse is given by 
P x = bT and P y = —aT (with T = —2ttk\) and clearly 
provides a measure of position through the guiding-center 
coordinates R x = —P y /T = a and R y — P x /T = b. 
Similarly, the angular momentum L of a static vortex is 
nonzero because the topological vorticity 7 is of definite 
sign. L is rather a measure of vortex size. For a single 
isolated vortex located at the origin we find numerically 
Ly = — 8.5kA. 

The linear momentum (impulse) of the AB pair shown 
in the first panel of Fig. Q] can be calculated approxi- 
mately if we suppose that the vorticity of each vortex is 
concentrated at its center. Considering that the antivor- 
tex has vorticity opposite to that of the vortex, we find 
Pab ~ {0,2nd), where d ~ 1/v is the relative distance 
between the vortex and the antivortex. A more precise 
calculation for v = 0.1 yields Pab — 62 

The single vortex C shown in the final panel of Fig. [T] is 
static and sits in the fourth quadrant at position (10, —6). 
This corresponds to an impulse P = (2ir x 6, 2ir x 10) ac- 
cording to Eqs. ©. The outgoing AB pair propagates 



with velocity v = 0.15 at an angle ip = 0.647T with re- 
spect to its initial direction of propagation. The impulse 
of a pair with v = 0.15 can be found in the appropri- 
ate table and figure of Ref. [13| to be Pab = 40, or, 
given the direction of propagation, Pab = (—36,-17). 
The result is that the AB pair changes its linear mo- 
mentum by APab = (—36, —56) due to scattering. The 
corresponding quantity for the single vortex is AP^ = 
(2ir x 6, 2tt x 10) = (38, 63). Conservation of linear mo- 
mentum would require that 

AP AB + AP C = (11) 

in the case of elastic scattering. This relation is approxi- 
mately satisfied in our simulation - the difference is due 
to spin-wave radiation produced during scattering. 

The numerical simulation of Fig. [T] shows that the lin- 
ear momentum (impulse) of a VA pair is transferred to 
vortex impulse which is tantamount to a translation in 
the vortex position. A genuine transmutation of momen- 
tum into position takes place during the three- vortex col- 
lision. Therefore, the definition of the impulse given in 
Eq. (J9j) is not only consistent with the symplectic struc- 
ture of the Landau-Lifshitz equation [lj| but it is indeed 
physically relevant. 

We may now use as the initial VA pair a solitary wave 
with larger velocity. We repeat the numerical simulation 
using a VA pair with v = 0.2. The time evolution is 
similar to that of Fig. [T] The velocity of the outgoing 
VA pair is v — 0.46 and the pair is scattered at an angle 
ip = 0.607T relative to the initial direction of propagation. 
The velocity of the outgoing VA pair is significantly larger 
that that of the incoming pair and indicates a significant 
inelastic component of the process. The position of the 
single vortex C after scattering is (5, —2). One can check 
that the difference in the VA pair linear momentum is 
transferred to the single vortex; i.e., to its translated po- 
sition. 

A simple approach to the dynamics of three vortices 
and specifically to the phenomena described in this sec- 
tion can be obtained through a collective collective- 
coordinate approach which is valid for sufficiently low 
velocity of the incoming VA pair. This approach is de- 
scribed in detail in the Appendix. In particular, we are 
able to deduce theoretically an approximation for the 
scattering angle of the VA pair. 

IV. SWITCHING OF VORTEX POLARITY 

The preceding numerical experiment was repeated for 
a Kelvin pair with relatively large velocity v = 0.5 for 
which the vortex and the antivortex are tightly bound 
at a relative distance d = 2.6 [HI]. The process is again 
illustrated by three characteristic snapshots in Figure [2] 
While the initial stages of the process are similar to those 
encountered in the case of slow Kelvin motion (Figure [T]) 
a substantial departure occurs when the pair now ap- 
proaches the target vortex. In particular, as soon as an- 



FIG. 2: Same as Fig. [T]but for a larger initial velocity v = 
0.5 of the AB pair. During collision, antivortex B begins to 
rotate around vortex C but the rotating BC pair is eventually 
annihilated leaving behind vortex A (with polarity opposite 
to that of the target vortex C) and a burst of spin waves that 
propagate away from the scattering region. 



tivortex B=(— 1,1) begins to rotate around the target 
vortex C=(l, — 1) they collide and undergo a spectacu- 
lar AM = 1 transition (annihilation) leaving behind the 
A=(l, 1) vortex which may be thought of as the target 
vortex C with polarity flipped from — 1 to 1 (vortex core 
switching) and a burst of spin waves propagating away 
from the scattering region. 

A crucial element in the preceding numerical experi- 
ment is the transient formation of a rotating VA pair (the 
BC pair) which is characterized by a nonzero skyrmion 
number J\f = 1. Nevertheless, the BC pair is annihilated 
during collision and thus leads to a topologically forbid- 
den transition that changes the skyrmion number by one 
unit or, equivalently, the polarity of the original vortex is 
flipped from —1 to +1. Such a process is indeed possible 
because a topologically nontrivial (and thus rotating) VA 
pair may shrink without encountering an energy barrier 
and eventually be annihilated when its size approaches 



FIG. 3: Same as Fig. [T] but for a solitary wave with ini- 
tial velocity v — 0.9. During collision the solitary wave is 
transformed into a VA pair. After collision the initial vortex 
has been annihilated with the generated antivortex while the 
generated vortex remains static near the origin. The net re- 
sult may again be thought as the original (target) vortex with 
polarity flipped from —1 to +1. 



the spacing of the underlying lattice. An excess amount 
of energy equal to 4tt is then released into the system in 
the form of spin waves [l4j |. 

According to the results in Ref. [HI solitary waves with 
velocities 0.78 < v < 1 apparently lose their vortex- 
antivortex character. Rather they are droplets of nonuni- 
form magnetization. In particular, the magnetization 
does not acquire the values m = (0,0, ±1) at any point, 
thus no vortex centers can be identified. In the remainder 
of this section we shall present simulations of the interac- 
tion of such solitary waves with a target vortex. We first 
choose an incoming solitary wave with velocity v = 0.9 
which may indeed be viewed as a weak disturbance of 
the ferromagnetic vacuum, as shown in the first panel 
of Fig. [3] As the solitary wave collides with the target 
vortex at the origin, it is transformed into a vortex and 
an antivortex, as seen in the second panel of Fig. [3J The 



FIG. 4: The second snapshot of Fig.[3]magnified. The VA pair 
created out of the incoming solitary wave is clearly shown 
here. The generated antivortex is about to be annihilated 
together with the original (target) vortex. 

same snapshot is shown magnified in Fig. |4j The antivor- 
tex is then annihilated with the target vortex and finally 
a vortex remains with polarity opposite to the initial one. 

Finally, we present a simulation of the scattering of a 
very fast solitary wave with velocity v = 0.95 which is 
close to the magnon velocity v = 1. Again, this wave 
does not display an apparent vortex-antivortex charac- 
ter, as can be seen in the first panel of Fig. [5] The soli- 
tary wave appears to shrink as it approaches the vortex 
at the origin and is deflected into the second quadrant. 
The deflected solitary wave is apparently not identical to 
the initial one due to energy radiation during scattering. 
The small deflection angle indicates a small change of the 
linear momentum during scattering, even though this is 
difficult to quantify accurately from the numerical data. 
But the observed small translation of the vortex from the 
origin after scattering indicates a correspondingly small 
exchange of impulse between the solitary wave and the 
vortex. No polarity switching occurs in this case. 

To summarize, a detailed numerical investigation of 
the three-vortex process for Kelvin waves with velocities 
in the allowed range < v < 1 suggests the existence of 
three characteristic regions separated by two critical ve- 
locities vi — 0.3 and — 0.9 such that < V\ < V2 < 1. 
For < v < V\, the Kelvin pair undergoes nearly elas- 
tic scattering of the type depicted in Figure [TJ For 
«i < v < V2, the process leads to a topologically for- 
bidden AAf — 1 transition of the type illustrated in Fig- 
ure [2] and Figure [3] There is also some evidence that 
fast Kelvin waves with velocities in the narrow range 
i>2 < v < 1 undergo a nearly elastic scattering without in- 
version of the polarity of the target vortex, as illustrated 
in Figure [5l 



V. CONCLUDING REMARKS 

A direct link between topology and dynamics has been 
known to exist in ferromagnetic media. In particular, 
the conservation laws of linear and angular momenta are 
profoundly affected by the underlying topology and ex- 
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FIG. 5: Three snapshots for the collision of a solitary wave 
initially located at (0, —15) and propagating with velocity 
v — 0.95, against a target vortex initially located at the ori- 
gin. During collision, the solitary wave as well as the vortex 
shrink, and finally the solitary wave is deflected in the sec- 
ond quadrant. The target vortex remains unaffected except 
that it is shifted slightly to a new position in the first quad- 
rant thanks to transmutation of solitary wave momentum into 
vortex position. 



plain some of the unusual d yna mical features of topo- 
logical magnetic solitons fl9l. I20I]. For example, a single 
magnetic vortex cannot move freely but is always spon- 
taneously pinned around a fixed guiding center. In con- 
trast, VA pairs with zero skyrmion number (M = 0) can 
move rigidly with constant velocity, in analogy with the 
Kelvin motion encountered in fluid dynamics jl3j . Free 
translational motion is again inhibited for VA pairs with 
nonzero skyrmion number J\f. Instead, such pairs un- 
dergo rotational motion around a fixed guiding center 

In the present paper we have taken the argument much 
further. Specifically, we have studied a three-body colli- 
sion during which a VA pair in Kelvin motion is scattered 
off an isolated single vortex. The analysis provided in- 
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teresting results at low velocities of the incoming Kelvin 
pair where the process is semi-elastic. The total skyrmion 
number is then preserved and the outcome of the collision 
is again a Kelvin pair and an isolated vortex. But the 
overall configuration after collision appears to be rather 
strange from the point of view of ordinary scattering the- 
ory. Specifically, the outgoing pair is scattered at a pe- 
culiar angle while the isolated vortex comes to rest at 
a different location. These facts would appear to vio- 
late conservation of linear momentum. However, the ob- 
served transmutation of momentum into position is fully 
consistent with the conservation laws quoted in Sec. IIIII 
and is indeed due to their special topological structure. 
Also note that an elegant mechanical analog of the three- 
vortex collision is provided by the analytical solution in 
terms of collective coordinates given in the Appendix. 
Finally, it is now evident that momentum-position trans- 
mutation is a generic and physically relevant feature of 
multivortex collisions and is thus the main new result of 
this paper. 

The process of three- vortex collision is also thought to 
be responsible for the vortex polarity switching observed 
in ferromagnetic elements [Til. 1T3| which entails a 
change of the skyrmion number by one unit. A mech- 
anism that changes the topological number of a mag- 
netic configuration makes it possible to obtain a con- 
trolled switching between topologically distinct (and thus 
robust) magnetic states, an issue of obvious interest for 
practical applications. Indeed, our numerical simulations 
within the 2D Landau-Lifshitz equation confirm the pos- 
sibility of polarity switching when an incoming VA pair 
in Kelvin motion with sufficiently high velocity collides 
with an isolated vortex. The process is highly inelastic 
in that a transient rotating VA pair is formed and subse- 
quently annihilated, leading to a change of the skyrmion 
number by one unit. This unusual phenomenon is due to 
the fact that a topologically nontrivial (and thus rotat- 
ing) VA pair can shrink during collision without encoun- 
tering an energy barrier and be annihilated when its size 
approaches the lattice spacing [l4[ ■ 



motion (HUH 



, * ,'dRi - 

I —TT X Z 



Ri — Rj 
\Rt - RA 2 ' 



(Al) 



We note that the mutual forces on the rhs depend only 
on the vortex numbers (re) and are attractive for pairs 
AB and BC and repulsive for AC. In contrast, the lhs of 
Eq. (|Alj) depends on the vortex skyrmion number and 
thus on both the vortex number (re) and the polarity (A). 
A more explicit form of Eq. (|A1[) is given by 



A, 



dt 



dXj 
dt 



\ Ri — Rj 



iRi-Rjl 2 ' 



(A2) 



applied for i=l,2 and 3, which leads to a system of six 
first-order differential equations that must be solved with 
initial condition 



X l = -d/2, X 2 = d/2, A 3 = 
Y x = -h , Y 2 = -h , y 3 = o 



(A3) 



which corresponds to the scattering process studied in 
Sec. IIIII In words, the topologically trivial (Af = 0) AB 
pair begins to move from a distance ho along the negative 
y axis with velocity v = l/d against a vortex C initially 
located at the origin. The initial configuration of the 
three vortices A, B and C is shown in the first panel 
of Fig. [6] which makes it clear that the triangle ABC is 
initially isosceles (AC=BC). 

The equations of motion (|A1[) and (|A2[) lead to four 
conservation laws; namely, the energy 



E = 2tt In 



\Ri — Rj] \R2 — Ri 
\Ri-R 3 \ 



(A4) 



the two components of linear momentum (impulse) 

P, = 2tt (Yi - Y 2 - Y 3 ) 

P y = -2ir{X 1 -X 2 -X 3 ). (A5) 
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and the angular momentum 

L = tt [-{Xl + Y 2 ) + (Xl + Y 2 ) + (A 3 2 + V 3 2 )] . (A6) 

The actual values of the preceding conserved quantities 
may be calculated from the initial configuration defined 
in Eq. (|A3|) . namely 



APPENDIX A: COLLECTIVE COORDINATES 

The main results of Sec. IIIII can be understood within 
a formulation in terms of collective coordinates where the 
three vortices A=(l, 1), B=(— 1,1) and C=(l, — 1) are 
approximated as pointlike particles whose coordinates 
Ri = (Xj,Yj), with i — 1,2,3, satisfy the equations of 



E = 2nhi(d), P x = 0, 



P y = 2nd, 



L = 0. 



(A7) 



However the above four conservation laws are not suffi- 
cient to readily solve the six differential equations con- 
tained in Eq. (|A"2|) . 

Hence, we first solved the equations of motion (|A2|) nu- 
merically with initial condition given by Eq. (|A3j) . The 
main features of the solution can be summarized as fol- 
lows: 
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FIG. 6: Four characteristic snapshots of the three- vortex colli- 
sion calculated within a collective-coordinate approach start- 
ing from an initial configuration with ho — 50 and d = 10. 
The three (anti)vortices form an isosceles triangle (AC=BC) 
with base of constant length AB=d at all times. See the Ap- 
pendix for further explanations. 



1. The triangle ABC remains isosceles (AC=BC) at all 
times, while the length of the base AB remains constant. 
Thus we may write 



AB = d, 



AC = DC = \lh? + ^ 



(A8) 



where only the height of the triangle depends on time, 
h = h(t), and satisfies the initial condition h[t — 0) = ho- 

2. There is a characteristic time to at which h(t = t ) = 
and the triangle degenerates into a straight line element 
ACB with vortex C located at a common distance d/2 
from both A and B. Therefore, to is the instance of closest 
encounter of all three (anti)vortices. 

3. For t > t , C emerges from the other side and 
forms again an isosceles triangle ABC. In the meantime, 
antivortex B encircles vortex C while keeping a constant 
distance d from its Kelvin partner A. Eventually, a new 
Kelvin pair AB is formed which moves off to infinity at 
a scattering angle in the third quadrant with asymptotic 
velocity equal to v = 1/d. 

4. Throughout the process vortex C moves on the lower 
semicircle of the circle 



(X 3 -d) 2 +Y 2 



(A9) 



In fact, it covers only a portion of the lower semicircle, 
up to a maximum angle that depends on ho and d. 
The preceding observations are consistent with the 



conservation laws and suggest the following parametriza- 
tion of the six coordinates: 

Xi = d — ^ cos ip + h sin ip, Y\ = — — sin^> — hcosip 



X2 = d — — cosC + hsintp, Y2 = - 
X 3 = d (1 — cos ip), Y3 = —d sintp 



— sin ip — h cos ip 
(A10) 



in terms of only two variables; namely, the height of the 
triangle h = h(t) and the scattering angle ip = ip(t) which 
is the angle between the symmetry axis of the triangle 
and the y axis. Then all six equations in (|A2[) reduce to 
the pair of equations 



A a = h 2 



~d L 

i. 

4 ' 



(All) 



where the overdot denotes time derivative, which must 
be solved with initial condition h(t = 0) = ho and tp(t = 
0) = 0. 

The second equation in (|A1 1[) can be integrated to 
yield 



1 



t 

2 ' 



to ~ id 1 
to = id 2 



h 

yd 

ho_ 

yd 



arctan 



arctan 



h 

yd 

ho_ 

yd 



(A12) 



Then we combine both equations in (|A1 1|1 to write 
dh/dijj = — (l/tf)(A 2 + d 2 ) which can also be integrated 
to yield 

— = tan[7(^ - VOL -A =tan[7V»o]. (A13) 
yd yd 

In words, at t = 0, h = ho and ip = 0. At t = to, h van- 
ishes and tj) reaches the value tpo ■ For t > to , h reemerges 
on the other side and assumes negative values in our con- 
ventions. In the far future (t — > 00), h ~ — t/d — > — 00 
and the scattering angle reaches the maximum value de- 
fined from 7(^0 - "0max) = —it/2 or 



In particular, if ho 



Ipmax = Co + 7T- 
2y 

00, Co ~ * 7r /(27) and 



7T 2vr 

— = — =, (as h 

1 V5 



00 



(A14) 



(A15) 



The preceding results are partially illustrated in Fig. [6] 
for an initial configuration specified by ho = 50, d = 10 
and ip = 0. 

In order to compare with the results on vortex scatter- 
ing presented in Fig. Q]we choose ho = 15, d = 9.5 and 
calculate t — 65.4, ipo = 0.854 and ^max = 0.727T. In 
the case of the simulation shown in Fig. [1] we have found 
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V'max = 0.647T. This number as well as the overall pic- 
ture are in good agreement with the results obtained by 
the collective coordinate approach. The observed small 
differences should be mainly due to the radiation present 
in the actual simulations and are expected to increase 



at larger velocity of the incoming Kelvin (AB) pair. In- 
deed, while agreement with the collective-coordinate ap- 
proach persists for u=0.2, the picture changes drastically 
for higher velocities, as discussed in Sec. IIVI 
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